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Abstract Universality and crossover is described for attractive and repulsive inter- 
actions where, respectively, the BCS-BEC crossover takes place and a ferromag- 
netic phase transition is claimed. Crossovers are also described for optical lattices 
and multicomponent systems. The crossovers, universal parameters and phase tran- 
sitions are described within the Leggett and NSR models and calculated in detail 
within the Jastrow-Slater approximation. The physics of ultracold Fermi atoms is 
applied to neutron, nuclear and quark matter, nuclei and electrons in solids when- 
ever possible. Specifically, the differences between optical lattices and cuprates is 
discussed w.r.t. antiferromagnetic, d-wave superfiuid phases and phase separation. 



1 Introduction and the Bertsch problem 

A decade ago at the 10th Manybody Conference G. Bertsch posed the problem (see 
IB for full text): 

How does a system of Fermi particles with infinite s-wave scattering length 
but vanishing interaction range behave? 



This seemingly innocent question turned out to be rather fundamental and trig- 
gered an explosion of interest and number of papers on physics now referred to as 
universal physics, the unitarity limit, BCS-BEC crossover, and strongly interacting 
and correlated systems in general. His question was partly motivated from nuclear 
physics and the physics of neutron stars where dilute gases of neutrons exist in the 
inner crust. The neutron-neutron 1 Sq scattering length a ~ — 18.5/m is long (and 
negative) compared to the order of magnitude smaller interaction range R ~ lfm, 
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because two neutrons almost have a bound state. Admittably, Bertsch was also mo- 
tivated by expectations that cold gases of Fermi atoms might be created like the 
BEC a few years earlier, and that it might be possible to tune interactions near Fesh- 
bach resonances between two hyperfine states in order to make the scattering length 
truly go to ±°o. Both were rapidly and successfully accomplished in a number of 
remarkable experiments where also the BCS-BEC crossover, multicomponent sys- 
tems, optical lattices, etc. have been studied. 

The solution to the Bertsch problem for two-component systems is as simple as 
it is fundamental |2][3): 

Since R — > and a — > ±°° both parameters must vanish from the problem leav- 
ing only one remaining length scale: the interparticle distance or equivalently 
the inverse of the Fermi wavenumber kp 1 . All thermodynamic quantities be- 
come universal. 

The proof is a dimensional argument which will be described in detail below. As 
an example, the energy per particle is an universal constant times the energy of a 
non-interacting Fermi gas (3/5)£> = 3h 2 k F /10m. The a — > ±°° limit is referred to 
as the universal or unitarity limit and the gas as a unitary gas. The name comes from 
scattering theory where e.g. the s-wave scattering cross section at relative momen- 
tum k is limited by Co <4-n/k 2 due to unitarity of the scattering matrix. For weak 
interactions Ob = 4-na 2 which would diverge like the energy per particle (see Eq. 
(0) - if extrapolated to a — >• ±00. In the unitary gas the cross section and the ther- 
modynamic quantities are instead limited by unitarity and universality respectively. 

In the neutron gas \a\ is long but fixed and the unitarity limit is defined as 
\a\ ^> ro R, where ro = (3«/47r)~ 1//3 ~ kp 1 is the interparticle distance, kF the 
Fermi wavenumber at density n = k F /3n 2 for two components/spin states. Thus the 
unitarity limit ranges from nuclear saturation density hum — 3/47Tfl 3 ~0.15firr 3 , 
and more than three orders of magnitude down in density to ~ |a|~ 3 . In this region, 
e.g. the energy per particle scales as an universal constant times Ef. The unitarity 
limit became easily accessible experimentally as it became possible to tune atomic 
interaction strengths near Feshbach resonances Q, where a —¥ ±°° corresponding 
to two-atom bound states at threshold. The scattering length could also be extended 
to positive values corresponding to bound molecular states and for strongly bound 
states a — > 0+ even a molecular BEC. 

The Bertsch problem described above was originally intended for two compo- 
nent system as spin 1/2 neutrons or ultracold atoms in two hyperfine states but it 
is also relevant for a nucleon gas, since the neutron and a proton have the weakly 
bound state of deuterium and therefore a large positive scattering length a ~ +20 fm. 
Isospin symmetric nuclear matter, however, has two spin and two isospin states i.e. 
four components and does not have a unitarity limit. It is in fact unstable at subnu- 
clear matter densities where a liquid gas phase transition occurs. This was success- 
fully described in O where estimates showed that as Pauli blocking between same 
spins is effectively reduced as the number of spin states increases, the Pauli pres- 
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sure can no longer overcome the unitary attraction. Only two component systems 
are stable and provides a unitary gas. Three component systems are marginal and 
four and higher component systems are unstable. Such multicomponent systems are 
now studied, e.g. 6 Li with three hyperfine states J5), 137 Yb with six nuclear spin 
states @, and heteronuclear mixtures of 40 K and 6 Li Q. Multi-component systems 
have intriguing similarities with neutron, nuclear and quark matter. In the latter color 
superconductivity between the 2 spin, 8 color and 2-3 flavor states may occur (8). 

When the interactions are varied, e.g. near a Feshbach resonance, the thermo- 
dynamic quantities depend on the dimensionless parameter akp. All the thermo- 
dynamic quantities become universal functions of the crossover parameter akp or 
equivalently x = 1/ akp in the sense that they do not depend on the system whether 
a gas of atoms, neutrons, or any other Fermi particle as long as R <C |a|,ro. The uni- 
versality argument is intimately connected to the smooth approach to the unitarity 
limit and crossover. At finite temperatures the universal thermodynamic functions 
depend on T /Ep as well. A decade earlier pairing models |]9] [TOl [TT] had already 
described the crossover from BCS to the BEC limit and calculated pairing gaps, 
transition temperatures and chemical potentials. 

In these lecture notes universality, crossover and correlations will be described. 
Not only the BCS-BEC crossover in uniform system but also the repulsive "ferro- 
magnetic" crossover, in multicomponent systems, traps and lattices. Mostly Fermi 
atoms are discussed but applications to neutron, nuclear and quark matter, nuclei 
and electrons in solids are made wherever possible. 



2 Universality and Crossover 

Universality is usually understood as the universal scaling in the unitarity (Bertsch) 
limit where the interaction range (R — > 0) and scattering length (a — > ±°°) are suffi- 
ciently short and large respectively that the thermodynamic quantities are indepen- 
dent of either or any other details of the interaction (such as effective ranges, etc.) 
and only given by universal constants. 

When only the interaction range is negligible (R — > 0) but the scattering length 
varies, as it is possible in ultracold atomic systems around Feshbach resonances, 
there is another dimensionless parameter in the problem usually taken as akp or x = 
1 1 akp. The thermodynamic quantities are thus universal functions of this crossover 
parameter varying from weakly interacting x — > — °° over the unitarity limit x = 
to the strongly interacting molecular limit x — > +°°. In pairing models IfTUl ITTI 
these limits are named the BCS and BEC limits respectively and the crossover as 
BCS-BEC crossover. When finite temperatures are considered the thermodynamic 
quantities are universal functions of both x and T/Ep lfT2l . 

It should be emphasized that there are several crossovers (see Fig. (Q~|i. The most 
common one is that for increasingly attractive scattering lengths also known from 
pairing models flOl fTTI and referred to as the BCS-BEC crossover. Here the cold 
atoms start out in the weakly attractive limit with BCS pairing i.e. x = 1 /akp — > — °° 



H. Heiselberg 



BCS (a-» 0_ ) Unitarity (a-> - 



mBEC 





o 


• 











• 


• 






o 




• 


o 





• 


o 









• 




o 


• 


• 


o 

















• 


•• 


o 





Dilute (a-> 



Gutzwiller 



Ferromagnet (a-> +■> 








• 




o 






• 


• 






o 




• 


o 





o 


o 




• 




• 


o 








o 




• 




• 





• 


o 




o 






• 


o 






• 




• 


o 





Fig. 1 Sketch of correlations in a two-component gas. Top row shows the gradual pairing of differ- 
ent components (spin states) with increasing attraction, i.e. the crossover from BCS via unitary gas 
to the molecular BEC. Bottow row shows the crossover for repulsive interactions from the dilute 
gas to either a hard sphere (Gutzwiller) or a ferromagnetic phase (see text). 



and goes over the unitarity limit x = to the strongly attractive limit x +°°, where 
a molecular BEC forms. Another crossover starts out with weak repulsive interac- 
tions, i.e. a small positive scattering length, and increasing it towards akf +°°. A 
recent experiment claims to observe a phase transition in a 6 Li gas to a ferromag- 
netic state for akF^l |fl3l . These two crossover are different because the latter has 
a node in the short range correlation function l22l . In principle there are a novel 
crossover for each number of nodes, however, they become increasingly short lived 
due to three-body losses. The different crossovers have different universal thermo- 
dynamic functions, universal parameters and quite different superfluid and ferro- 
magnetic phases. 

We will in the following subsections describe universality and crossover mainly 
for the BCS-BEC crossover and briefly for the "ferromagnetic" crossover. A number 
of thermodynamic functions are described which are calculated within the Jastrow- 
Slater approximation for the many-body wave function. 
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2.1 Universality in the unitarity limit 



The dimensional argument implies that any thermodynamic function, for example 
the energy per particle in units of the Fermi energy, can at zero temperature only 
depend on one dimensionless scale such as x 

| = |ifr(l +£(*)). (1) 

Here, the ratio between the interaction and kinetic energy j3(x) is a "universal" 
many-body function in sense that it and all other thermodynamic functions are in- 
dependent of the system, e.g. dilute neutron matter, helium gases, atomic gases, etc. 
Earlier definitions are tj = (3/5)(l +J3) |14) and v = (5/3)0 0. 

The universality and crossover hypothesis implicitly assumes that the transition 
to the unitarity limit and the crossover is smooth, i.e. the derivatives of fi(x) are 
finite. Taking the derivative with respect to density gives 

^ = 5(i+^W)^-^'W (2) 

an 5 an 5n 

where j3'(x) = df5(x)/dx. We observe that derivatives of j3 w.r.t. densities or : 
kpd /dkf = —xd/dx, always brings a factor x or 1 /a. In the unitarity limit x = 0, 
such terms as the second term in Eq. (O therefore vanishes and only the universal 
parameter /3(0) remains. By repeating the argument all density derivatives of the 
energy such as chemical potential, pressure, first sound speed, compressibility, etc. 
depend only on this one universal parameter j3 (0). Other thermodynamic variables 
than density derivatives of E/N can, however, depend on higher derivatives /3 W (0) 
as we shall see below. 

The first estimates based on Pade' approximants and Galitskii resummation gave 
J3(0) ~ —0.67 [2, 3 1 in the unitary limit whereas the Jastrow-Slater approximation 
yielded j3 (0) ~ —0.54. The Leggett pairing model, which will be described in detail 
in a later section, predicts j3 (0) ~ —0.40. Ground state energies have been calcu- 
lated more accurately numerically by Monte Carlo for systems with a finite num- 
ber of Fermions in a box. The earliest actually use the Jastrow-Slater as trial wave 
functions for Green's function lfl4l and diffusion (T5 1 Monte Carlo. Since pairing 
becomes important at crossover a more general BCS wave function based on the 
Jastrow wave function was also employed. Extrapolating to a large number of par- 
ticles they obtained J3 (0) = -0.56 ± 0.01 and J3 (0) = -0.58 ± 0.01 respectively in 
the unitary limit. 

Recent measurements J4] [16] [T7] [18] [19) of j3(0) confirm the unitarity limit near 
Feshbach resonances. Several experiments with trapped Fermi atoms have recently 
measured energies in the strongly interacting or dense limit near Feshbach reso- 
nances. The energy in the trap (excluding that from the harmonic oscillator poten- 
tial) is E/N = (3/8)£>0~+/K°l where e f = (3A0 1/3 /iO) is the Fermi energy in 
a trapped non-interacting gas. The first measurements by the Duke group 0] mea- 
sured the energy of 6 Li Fermi atoms near a Feshbach resonance from expansion 
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energies. These early measurements were later corrected for thermal energies and 
find JS (0) = -0.4 ±0.1. HE) With the discovery of a molecular BEC the Inns- 
bruck group lfl9l has been able to measure the size of the atomic cloud, which scales 
with (1 + /3) 1 / 4 , around the Feshbach resonance at very low temperatures, and find 
/3 (0) = -0.68 ± 0. 1 . Recent accurate measurements find J3 (0) = -0.49 ± 0.02. ED 

Other thermodynamic variables may in the unitarity limit depend on other uni- 
versal parameters. For example, the number of closed channel molecules is propor- 
tional to P'(x) [32 33) and one obtains J3'(0) ~ —0.5. From the slope of collective 
modes near the unitarity limit J3'(0) ~ — 1.0 J4][l7). Generally, the universal func- 
tion j3 (x) is given by a Taylor expansion in terms of an infinite number number of 
derivatives j3 W (0), which all are universal parameters that determine all the univer- 
sal thermodynamic functions. 

Unlike neutrons, ultracold atoms can have many internal bound states and cor- 
responding Feshbach resonances as the background magnetic field is increased. 
Thus solutions are multivalued for a given scattering length where the gas is in a 
metastable state. For example, starting from a non-interacting gas (a = 0) and in- 
creasing the scattering length a > we also approach an unitarity limit a —> +°° 
l2ll . which differs from that in the BCS-BEC crossover when the scattering length 
is decreased a < 0. This unitarity limit is similar to that for bosons as a — >• +°° 11221 . 
where the scattering length must be positive in order for the system to be stable. 
As will be described in detail the two-body wave function has a node and therefore 
the universal functions differs from those in the BCS-BEC crossover. In principle 
a new universal limit and crossover exist for each number of nodes, n = 1,2,3, 
in the two-body correlation function. These will, however, be increasingly unstable 
towards three-body losses. The « = 1 has been observed for bosons in the atomic- 
molecular transition in a 85 Rb BEC 0231 and is in agreement with the predicted value 
from JS with K\ = 2.80. For fermions the predicted value is calculated below within 
the Jastrow-Slater approximations /3i(0) ~ 2.93 which is compatible with a recent 
experiment |fl3l . This value is about six times |J3(0)| demonstrating that these two 
universal limits are very different. 



2.2 Thermodynamic crossover functions 

As argued above any thermodynamic function can only depend on one dimension- 
less scale such as x at zero temperature when the range of interaction is sufficiently 
small. For example, the energy per particle in units of the Fermi energy is given in 
Eq. (Q]i in terms of the universal many -body crossover function j3 (x) = Emt/Ekin- It 
is well known in the dilute limit for Fermions l24l 



N 



E 




an + ... 



(3) 



to leading orders, whereas for molecular bosons 
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Fig. 2 The thermodynamic functions J3(*) and /3i(jc) within JS. Also shown are the BCS limit 
of Eq. $5$ and the chemical potential in the Leggett crossover model (see text). On the BEC side 
(x > 0) the molecular energy been subtracted. 
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where the molecular scattering length is ~ 0.6a |54l[l6l . Thus we obtain 

( 10/(9ttjc) +4(11 -21n2)/(2l7T 2 x 2 )- 

P(x) = lp(0)+P'(0>+- 

[ -5x 2 /3 - 1 + {a M /a)5/(\%%x) + .. 

In Fig. d2J various model calculations of the universal function j3 (x) are shown. 
The chemical potential is in terms of the universal function 

fdE 



(4) 




(5) 



v.s 



E P [l+p--xP' 



(6) 



In both the hydrodynamic limit and for a superfluid gas the first sound is given 
by the adiabatic sound speed 



m \ on 



1 



v,s 



(7) 



where Vf = Kkp / n 
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The compressibility k = n 2 (dn/dji) is related to first sound as 



K- _1 =n 2 (^| =mnct. (8) 
\ dn Jv,s 

The pressure P = n(fl —E/N) = n 2 dE/dV is 

P{x) = ^E F n[l+P-xP'/2]. (9) 

The polytropic index P °c is defined as the logarithmic derivative of the 
pressure 

Y(x) -Pdn l + P-xp>/2 ■ ( ' 

7 approaches 2/3 in both the dilute and unitarity limit whereas it approaches 7=1 
in the molecular limit. This index determines the frequency of collective modes in 
traps 1601 . 

We notice that all the above thermodynamic derivatives w.r.t. density only de- 
pend on one universal parameter j3(0) in the unitarity limit x = as argued above. 
Most thermodynamic quantities behave in the unitary gas as in a free Fermi gas ex- 
cept for the universal scaling factor (1 +J3(0)). This factor, however, cancels in the 
polytropic index and therefore 7=2/3 in both the unitarity and BCS limit and the 
collective modes are the same as has been verified experimentally J4). 



2.3 Finite temperature 

At finite temperature the thermodynamic functions also depend on the parameter 
T/E F (see, e.g., G] [12] |25] |26)). Usin 8 (d(E /N)/dT) v , N = T{ds/dT) v . N = cy, 
where s = S/N is the entropy and cy the specific heat per particle, we obtain from 
Eq. QJ 

P( x ,T/E F )=P(x,0) + 5 -j Q T c v (x 1 ^-)^-. (11) 

The entropy is at temperatures well below the superfluid transition temperature 
T <C T c given by phonon fluctuations 

,= *4( T -)\ (.2) 
45« \c s ) 



Inserting in eq. ( fTTT i gives at low temperatures. 
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P(x,T/E F )=P(x,0) + 5^-s. (13) 

The temperature corrections scales as ~ (T/Tpy and since T^jT c ~ 0.19£>, the 
temperature dependence is almost flat as also observed in |20l . 

In a Fermi liquid s = (n m* j2m)T JEf, where m* is the effective mass, and thus 

^(jc > r/£F)~j5(x,r c /^) + ^ -g-s (14) 

when r c < r <C £V . Recent accurate measurements find j3(0,0) = —0.49 ±0.01 
and m*/m— 1.13 ±0.03 ||20l in the unitarity limit. These allow us to extract the 
Landau parameters at zero temperature = (1 + J3(0,0))m*/m — 1 ~ —0.42 and 
F( = 3(m*/m- 1) ~ 0.39. In the BCS limit m*/m = 1 + [8(71n2 - 1)/I5n 2 }a 2 kj- 
andF = l + (10/9s;)fljt F . 

As the temperature decreases towards the critical temperature the change in /3 (T) 
from a Fermi liquid Eq. ( fT~4T > towards a superfluid Eq. ( fT~3T > allowed determination of 
TrlEr = 0.19 ± 0.02 ll20T 

At temperatures below T < T t the normal and superfluid components lead to two 
(first and second) sound modes in the collisional limit. In the BCS and BEC limit 
these are also referred to as the compressional and thermal sound modes. Their 
velocities u\ and ui are given by the positive and negative solutions respectively of 




u 2 = tS^2 ± ^tS^l) _ c l c 2. (15) 

The thermodynamic quantities entering are the adiabatic c| = (dP/dn)s and the 
isothermal c\ = (dP/dn)r compressional sound speed squared. The "thermal" 
sound wave c\ = n s s 2 T /n n cy acts as a coupling or mixing term and is small in 
the two limits. Here, n = n n + n s is the total, n„ the normal and n s the superfluid 
density. The difference between the adiabatic and isothermal sound speed squared 
can also be expressed as 

ds\ 2 n 2 T 

7r) ' (16) 

an J T cy 

The mixing of the compressional and thermal sound modes is particular interesting 
at crossover where they mix and couple strongly. They undergo avoided crossing, 
i.e. the compressional change smoothly into a thermal sound mode and visa versa 
around the unitarity limit lf25l . 
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2.4 Jastrow-Slater approximation 



The Jastrow and Jastrow-Slater (JS) approximation was among the earliest models 
applied to the unitarity limit and crossover J3|. It has the advantage that it provides 
analytical results for the universal parameters that are easy to understand and has 
proven to be quite accurate when compared to experiment. It gives an ansatz for the 
strongly correlated wavefunctions which is also the starting point as trial wavefunc- 
tions in Monte Carlo calculations lUSl , Finally, JS can be generalized to describe 
both bosons and fermions with any number of spin states as well as other universal 
limits. 

The Jastrow and the JS approximation methods were developed for strongly in- 
teracting and correlated Bose and Fermi fluids respectively such as 4 He, 3 //e and 
nuclear matter in |29l and has more recently been applied to kaon condensation 
ll30ll . cold atomic Fermi |2T| and Bose 12211 gases. As explained in these references 
the JS wave function 

Wj S (r 1 ,...,r N ) = €> s ]J<l>(r i -r f ), d7) 

hf 

incorporates essential two-body correlations in the Jastrow function (j)(r). The an- 
tisymmetric Slater wave function <P$ f° r f ree fermions <P$ insures that same spins 
are spatially anti-symmetric. The Jastrow wave function only applies to particles 
with different spins (indicated by the primes). The pair correlation function (j)(r) 
can be determined variationally by minimizing the expectation value of the energy, 
E/N = (f \H\W) I (f|f), which may be calculated by Monte Carlo methods that 
are fairly well approximated by including only two-body clusters. The basic idea 
of this method is that at short distances r <C ro the Jastrow function (r) obeys the 
Schrodinger equation for a pair of particles interacting through a potential U (r) 



m dr 2 



r<j>(r)=e M r<Kr), (18) 



where the eigenvalue energy of two atoms Em = 2Ej„ t /N. Many-body effects be- 
come important, when r is comparable to ro, but are found to be small in lowest 
order constrained variation (LOCV ||29l [3] [30j HD). Here the boundary conditions 
that <p(r > d) is constant and <p'(r = d) = are imposed at the healing distance 
d, which is determined self consistently from number conservation (see l2TI for 
details) 

V-l f d (t) 2 (r) j 

-n \ -^)-p-Anr 2 dr = 1 (19) 



v Jo (p 2 (d) 

Note that a given component only interact and correlate with the (v — 1) other 
components which explains the prefactor. In the dilute limit (r) ~ 1 and so d = 
(v/(v — l)) 1 / 3 ^. In the unitary limit a — > ±°° the healing length now approaches 
d = r (2v/(v - 1)3)'/ 3 = (37r/(v - l)) 1/3 /t F '. Generally d ~ r . The boundary 
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Fig. 3 K 2 vs. scattering length in the JS model in units of interparticle distance ro- 



condition at short distances is given by the scattering length (r^>)'/r<j) = —l/a at 
r = 0. 



2.4.1 Attractive crossover 



For negative scattering lengths or negative energies Em < the solution to the 
Schrodinger equation gives a wave-function r<j>(r) sinh[£(r — b)} for weak inter- 
actions which change to r<j>(r) <=c cosh[A:(r — b)]) for stronger interactions (see Fig. 
(|4|). The boundary conditions and number conservation determine the phase kb, the 
energy and the healing length d. For small scattering lengths b = a whereas b = 
in the unitarity limit. The interaction energy Ej„,/N = —h 2 k 2 /2m is given by |f3T| 

a K~ 1 tanh K - 1 
d 1-KtanhfC 



(20) 



with K = kd. In the BCS limit Eq. ( f20b gives the correct interaction energy per 
particle, Eq. (O. In contrast the negative energy solution to Eq. ( f2Qb reduces in 
the unitarity limit to fetanhK - = 1 with solution Kb = 1.1997... As the scattering 
length cross over from — °° to +°° the negative energy state is analytically continued 
towards the molecular bound state with E/N = —h 2 /2ma 2 as a — > +0. 

In addition to the interaction energy Ej„ t = K 2 /2md 2 as calculated above a kinetic 
energy (3/5)£j? appears due to the Slater ground state. The total energy becomes 
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e 3 n 2 K 2 

N = 5 Ef ~2^- (21) 
From the definition Eq. © we obtain the universal function 

»=-\{£i- <22> 

In the unitarity limit k=k = 1.1997 and dk F = (37r/(v - l)) 1 / 3 and the universal 
parameter is for a general number of spin states 



5 /v-1 



2/3 



^ (0) = -3^J * ^ 

Thus in a two-component system j3(0) = —0.54 which lies between Monte Carlo 
results J3 (0) ~ -0.56 E] and recent experimental data j3 (0) ~ -0.49 l20l . 

The universal function J3(x) is shown in Fig. (f2]i for the JS approximation. Also 
shown is the chemical potential within the Leggett crossover model which will be 
discussed below in connection with pairing. The JS model includes self energies 
and is therefore a better approximation on the BCS side (x < 0). On the BEC side 
(x > 0) both models approach the molecular energy (subtracted in Fig. (f2]i) but to 
next orders both model overestimates the energy: the JS model by the Slater energy 
and the Leggett model by overestimating the molecular scattering length om by a 
factor ~ 3. 

The dependence on the scattering length is given by Eq. (l20l . By taking the 
derivative on both sides w.r.t. x we obtain 

« »=-y^F"-°- 48 - <24) 

This JS prediction is somewhat lower than that of Monte Carlo: j3'(0) = — 1.0 ±0.1 
lfl4l[T5l . The slopes of the axial and longitudinal collective frequencies of trapped 
unitary Fermi gases are directly proportional to jS'(0) and measurements also give 
J3'(0) ~ — 1 .0 [17) . Analysis of the number of closed channel molecules indicates 
J3'(0) ~ -0.5 |33][32)- The Leggett pairing model J3'(0) ~ -1.0. 

The universal parameter j3'(0) will be related to short range correlations in sec. 
(El- 



2.4.2 Repulsive crossover 

Starting from a positive scattering length the and positive interaction energy the 
solution to the Schrodinger equation gives a wave-function r0(r) sin[£(r — b)] 
(see Fig. dU). The boundary conditions and number conservation gives l22l 
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Fig. 4 The JS wave function r<j)(r)/d<j>(d) in the dilute limit (if = 0), the attractive unitarity limit 
(Kq), towards the mBEC (K = 2), the repulsive unitarity limit with one (Ki) and two nodes (K2). 



q = K-'tanK-l 
c/ 1 + iftan ff 

As in the BCS limit Eq. d25l l gives the correct interaction energy per particle in 
the dilute limit, Eq. ©. In the unitarity limit a — > +°°, the positive energy solution 
reduces to KtmK= —1 with multiple solutions K\ =2.798386.., K2 =6. 1212.., etc., 
and asymptotically k„ = nK for integer «. Generally, re = 0, 1,2,.. is the number 
of nodes in the JS wavefunction and each determines a new universal limit with 
universal parameters depending on the number of nodes. The phase in the wave 
function is kb = n{n — 1/2) whenever the unitarity limit of n nodes is encountered. 
The resulting energy is 

E 3 h 2 K 2 

N = 5 Ef + ^ (26) 



and we obtain from Eq. ([T|i 



'w-Ks) ■ <27) 



Because Eq. d25l ) have a string of solutions for given a, k and /3 are multivalued 
functions for a given scattering length or x, which we distinguish by an index j3 n 
referring to the number of nodes n = 0, 1,2, .... The attractive case discussed above 
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is the n = case, i.e. universal function in the BCS-BEC crossover is also j3 = Po- 
ll = 1 is the first weakly repulsive crossover where a ferromagnetic phase transition 
may occur and is also shown in Figs. (f2]) and ([3]). 

The universal parameter in this repulsive unitarity limit is 

Pi (x = 0) = +| (^fj k\ ~ 2.93(v - 1 ) 2 / 3 . (28) 

The index denotes that this universal parameter is for the wave function with one 
node. It has recently been measured in a 6 Li gas in two spin states ff3l . The chem- 
ical potential in the optical trap almost doubles going from the non-interacting to 
the unitarity limit. As shown below it scales as jj, oc jSi (0) giving j3i (0) ~ 3 

compatible with JS. However, the experiments also indicate a phase transition in 
repulsive unitarity limit from a paramagnetic to a ferromagnetic phase. This phase 
transition is supported by the JS models as will be shown in sec. ( 12.7b . 

As for the attractive crossover we calculate within JS the derivative from Eq. $25[ 

« (0) = -TW^-'- 38 ' < 29 > 



2.5 Short range correlations 

The short range two-body correlations are connected to thermodynamic quantities of 
the system J34) . We shall here derive one important example at zero temperature and 
use the Jastrow-Slater results above to actually calculate a new universal number. We 
refer to 11351 l26l for finite temperature extensions. 

Following l26l we scale the two-body potential by a factor A and obtain from 
Feynmann-Hellmann 

= V { JX ] = I d 3 rU(r){ V l(r) Y l(0)Y2(0)Yi(r))x. (30) 
Here the correlation function at short ranges is 

(Yl(v)Yl(0)Y 2 (0)Y l (v)) l = (l)l(r)=c(^-^j , (31) 

where <px is the solution to the two-body Schrodinger equation ( TT8l but with po- 
tential XU (r). C is a normalization factor such that (1 /V) J (r) = (n /2) 2 , and 
will later be interpreted as the correlation strength. 

Taking the derivative of the Schrodinger equation w.r.t. X, multiplying by r<j>i, 
and integrating over r we obtain 

r C drU(r)(r^(r)) 2 = --^, (32) 
Jo tn aX 
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where r c is any range between R<c c < tq. Combining Eqs. ( f30l > and (l32l gives 

dE/V 4nC 



da- 1 



(33) 



From Eq. ([]]) we can relate the correlation strength to the derivative of the universal 
functions j3„(x) as 

c » = -^a:m, (34) 

where the node index n distinguishes the various universal functions as discussed 
above. In the unitarity limits we can determine C n [x = 0) by inserting j3q(0) = —0.54 
or f}[ (0) = — 1 .38 from JS Eqs. d24l i and d29l l, etc. C n (x) can also be obtained directly 
from the JS wave-function, normalization and Eq. (l3TT l. 

The Fourier transform of the wave-function at large momenta is determined by 
the short range correlations <f)(k) = f d 3 re' k ' r (j)(r) ~ \K\fCjk 1 . Therefore the mo- 
mentum distribution for a spin state has the tail 

n a (k) = </> 2 (fc) = (4x) 2 C/k\ (35) 



at large momenta as is well known for Fermi liquids [24]. The correlation strength 
is directly measured in the number of closed channel molecules l32l 

dE da' 1 1 
da- 1 dB jl M ' 



Nm = ^t— , (36) 



where /Zm is twice the magnetic moment of the atom. The scattering length de- 
pends on the applied magnetic field as a ~ at, g [l — AB/ (B — Bo)] near a Feshbach 
resonance at Bo. Experiments find for the BCS-BEC crossover 1321 /3q(0) ~ —0.5 
which agrees well with that of JS and is compatible with model estimates 



2.6 Instability in multicomponent systems 

The original Bersch problem was intended for two-component systems as neutron 
matter or spin-balanced ultracold atoms in two hyperfine states. Early extrapolations 
of the dilute energy of Eq. (0 to a — > — °° led to the belief that the unitarity limit was 
unstable towards collapse |[36l as a nuclear gas whereas in fact the two-component 
system is the only stable unitary gas. The predicted scaling for the neutron gas at 
subnuclear density as given by the unitary gas of Eq. (Q]i triggered the memory of 
Vijay Pandharipande ll37l . who had the impression that Bethe and Brueckner were 
aware of this scaling when they looked at low density neutron matter in neutron 
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star context back in the '50s l38l |29l Fl Calculations by Carlson et al. |[T4ll have later 
confirmed this scaling in neutron gases. 

Multicomponent systems are now studied, e.g. 6 Li with three hyperfine states 
0, 137 Yb with six nuclear spin states |6|, and heteronuclear mixtures of 40 K and 
6 Li (7). Such multi-component systems have intriguing similarities with neutron, 
nuclear and quark matter where color superconductivity between the 2 spin, 8 color 
and 2-3 flavor states may occur JH). 

For the gas to be stable towards collapse in the unitarity limit the energy must be 
positive, i.e., 1 +/3(0) > 0. In the JS approximation we obtain in the unitary limit 
from Eq. <|23]l j3 (0) = 0, -0.54, -0.85, -1 .12, for v = 1 ,2, 3,4, ... spin states re- 
spectively. Therefore the JS approximation predicts that up to v < 3 spin states are 
stable in the unitary limit whereas in the Galitskii approximation only v = 1 , 2 are 
stable [3 1. Pauli blocking is effectively reduced in many component system and can 
only stabilize one, two and perhaps three component Fermi systems in the unitary 
limit. The stability of two spin states towards collapse has been confirmed for a 6 Li 
and 40 K gases near Feshbach resonances. The marginal case V = 3 has been studied 
with 6 Li atoms, which have three spin states with broad and close lying Feshbach 
resonances, and although the loss rate of atoms is large near Feshbach and Efimov 
resonances jfQ the gas is sufficiently long lived for measurements and does not col- 
lapse. It has long been known that neutron star matter |39l with two spin states 
likewise has positive energy at all densities whereas for symmetric nuclear matter 
with two spin and two isospin states, i.e. v = 4, the energy per particle is nega- 
tive. Nuclear matter is therefore unstable towards collapse and subsequent implo- 
sion, spinodal decomposition and fragmentation at subnuclear densities |40| . Above 
nuclear saturation densities, kpR^l, short range repulsion stabilizes matter up to 
maximum masses of neutron stars ~ 2.2M©, where gravitation makes such heavy 
neutron stars unstable towards collapse (39). The conjecture is therefore that the 
V > 3 Fermi systems are unstable and non-universal in the unitarity limit a —> — °°. 
Bose atoms corresponds to v = °°, since no Pauli blocking applies, and are unstable 
for negative scattering lengths. 

The v = 3 system is relevant for several reasons. Traps with ultracold 6 Li atoms 
with three hyperfine states are sufficiently stable and long lived @ to be studied in 
detail. The three Feshbach resonances are, however, separate in magnetic field so 
that the unitarity limit is not simultaneous for the three components. The three body 
system has interesting Efimov states which are non-universal ll4Tll42ll . i.e. besides 
the scattering length the system depends on an additional potential parameter such as 
the effective range. This non-universality persists for three bosons confined in a trap 
l43l and non-universality is therefore expected for a gas with three state Fermi atoms 
as experiments for the 6 Li system also indicates (5). Therefore three component 
systems do not have an universal limit and are only marginally stable, i.e. stable 
towards collapse but suffer three-body losses. 

In spin polarized systems or systems with different densities of the spin compo- 
nents the stability conditions depends on the various component densities. Also the 

1 However, no written reference has been found in the works of Bethe, Brueckner or Pandharipande 
so far. 
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system may undergo phase separation into a more symmetric phase and an asym- 
metric phase as, e.g., nuclei and a neutron gas in the inner crust of neutron stars. 
Similarly, spin polarized atoms at low temperature in traps may for strong attractive 
interactions separate into a paired spin balanced phase in the centre with a mantle 
of excess spin atoms. Likewise phase separation of strongly repulsive fermions may 
separate into domains of ferromagnetic phases in the centre with a paramagnetic 
mantle around. 



2. 7 Repulsive interactions and itinerant Ferromagnetism 

For a small positive scattering length the Fermi gas is a paramagnet (PM). For 
stronger repulsion Stoner l44l predicted a phase transition to a ferromagnet (FM) 
which a recent experiment claim to have observed lfT3l . Stoner's argument was 
based on the dilute equation of state of Eq. (O which generally for a spin polarized 
two-component system of total density n = n \ + «| and polarization rj = (n i — tvf) jn 
is 



E/N = ^E F 



20 

(l+7 7 ) 5 / 3 + (l-77) 5/3 + ^(l+'?)(l-'?)^« 



(37) 



Expanding for small polarization gives an equation of the Ginzburg-Landau type 

E/NE F ~ 1 + 2-ak F + Ul- -ak F )ri 2 + 3- 4 ri 4 + 0(ri 6 ). (38) 
5 in in 

It predicts a second order phase transition at ak F = n/2 from a PM to a FM with 
polarization r\ = ±y/27(2ak F — l/n). Due to the small fourth order coefficient it 
quickly leads to a locally fully polarized system rj = ± 1 . 

Unfortunately the predicted transition occurs close to the unitarity limit where the 
dilute equation of state is not valid. Higher orders may be important as exemplified 
by including the next order correction of order a 2 . It changes the transition from 
second to first order (45] at low temperatures up to a tri-critical point at T c ~ 0.27>, 
where the transition becomes second order again. However, the dilute expansion 
remains invalid in the unitarity limit. 

The Jastrow-Slater approximation extends to the unitarity limit also for positive 
scattering lengths as discussed above. Number conservation of the various spin den- 
sities is automatically included in the healing length, see e.g. Eq. < TT~9T >. As result the 
energy is 



E 3 
- = — E F 
N 10 * 



( i + 77 ) 5/3 + ( i - j? ) 5/3 + ft ( i + n ) ( i - n ) 2/3 + ft ( i + n ) 2/3 ( 1 - n ) 

0(ri 6 ). (39) 
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Pi{x, 77) is now also a function of polarization but for simplicity we shall ignore 
the dependency this dependency as we expect it to be minor. Expanding for small 
polarization we find 

E/NE F ~ |(1 +J3i) + 1(1 - Ipfo 2 + 3-\l - j3077 4 . (40) 

Truncating to order rj 4 would (erroneously) predict a second order phase transition 
from a PM to a FM at P\(x) = 5/7 ~ 0.71. The fourth order term is even smaller 
than the dilute prediction of Eq. d38l l. and it is therefore necessary to include higher 
orders. By equating the energy of the unpolarized gas, ~ (1 + Pi) with that of a 
fully polarized gas, ~ 2 2 / 3 , we find a first order transition at j3i = 2 2 / 3 — 1 ~ 0.59, 
since this value is smaller than 5/7. In view of the approximations made in this 
model calculation and the proximity of the two Pi values for the first and second 
order transitions, we can not reliably determine whether the order of the PM-FM 
transition is first or second. However, since the repulsive interaction energy in the 
unitarity limit /3i(0) ~ fcf/(37l) 2 ' 3 ~ 2.93. is much larger than the critical value 
Pi ~ 0.59 we can safely conclude that the transition to a FM does take place at a 
value corresponding to akp ~ 0.85. 

In a recent experiment the transition is observed around ak F ~ 2.2 at tempera- 
tures T/Tp =0.12 and ak F ~ 4.2 atT/Tp— 0.22 1131 . Presumably, the critical value 
for ak F is smaller at zero temperature. Also the Fermi wavenumber k F is the central 
value for a non-interacting gas which is larger than the average value over the trap 
of the gas that is further expanded due to repulsive interactions. More experiments 
will determine j3j (x), a possible ferromagnetic phase and critical value for akp. 

When the number of atoms in the two spin states is balanced the ferromagnetic 
domains of 77 = ± 1 coexist. Their domain sizes may be to small to observe in present 
experiments lfl3l . The densities of the two components will, however, have interest- 
ing distributions for unbalanced two-component systems in traps, where the mi- 
nority component will be suppressed in the centre and both phase separation and 
ferromagnetism can occur. In three component systems, when there are more than 
one Feshbach resonance as in 6 Li, with Feshbach magnetic field such that two res- 
onances oi2 and an are large but «23 small, the atoms will separate between a FM 
phase of 1 and a mixed FM phase of 2+3 with different densities. 

It should be emphasized that for positive scattering lengths the wave function 
and thus the correlations function between fermions of unlike spin and bosons % = 
r(j) ~ sin(kr — b) has a node somewhere within the interparticle distance [0;ro] (see 
ll22ll and Fig. It does not vanish as r — > as does the wave function for a short 
range repulsive potential as in hard sphere scattering, where a ~ R. Therefore the 
Gutzwiller approximation discussed in [46] applies to hard sphere gases, strongly 
correlated nuclear fluids and liquid helium but not to the repulsive unitarity limit of 
ultracold gases (see Fig. ([T). 
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3 Pairing in uniform systems and the BCS-BEC Crossover 



Experiments on superfluid Fermi gases have recenty confirmed the BCS to BEC 
crossover models J9] [10] QT| [56] [57) for pairing that was developed after the in- 
vention of BCS theory. The development of BCS-BEC crossover will be described 
historically with increasing level of complexity. 



3.1 BCS limit 



Bardeen, Cooper and Schriffer (BCS) first wrote down the famous gap equation l47l 
for an attractive two-body interaction U(r) < 



1 , . 1 - 2f{E„,) 



2E„ 



(41) 



which can be elegantly derived via the Bogoliubov transformation ll48l . Here, /(e) = 

(exp(e — ll)/T) + is the Fermi distribution function, = (e^ — jj,) 2 + Aq 

the quasi-particle energy and = h 2 k 2 /2m the free particle energy. In the BCS 
limit number conservation insures that ji = Ef. t/(p,p') = t/(p — p') is the Fourier 
transform of U (r). In metals phonons provide a small residual attractive interaction 
t/(p', p) ~ — V„t, with a cutoff of order the Debye frequency cod- As result we obtain 
from Eq. < f4TT > the BCS s-wave gap Aq = 2h(ODe~ x l N ^ v i> h at zero temperature. 

In the years immediately after BCS was developed the gap equation was gen- 
eralized in terms of scattering lengths [49]. When the interaction U(r) is short 
range, its Fourier transform U(p — p') is long range in momentum. It is then con- 
venient to replace the interaction by its scattering matrix T = U + UGqT, where 
Go = 1/ (2£k — iS), is the vacuum propagator for two particles. At low momenta 
the scattering matrix is given by the s-wave scattering length T = Anh 1 a/m = U$. 
Eliminating U in the gap equation (HTl i gives 



l= U0y 

2V^ 



j__ l-2/(£ k ) 



(42) 



Note that the difference between the vacuum and in medium Green's functions au- 
tomatically cuts off the high momenta, which are now included in the scattering 
length. Solving this gap equation at zero temperature for a Fermi gas interacting 
through an attractive scattering length a < gave a pairing gap in the dilute limit, 

\a\k P <:l, HO 



A = -jE F exp 



2akp 



(43) 
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Gorko\0 pointed out that many-body effects (induced interactions) lead to the next 
order correction in the interaction 150 l5"Tll52ll53l 

uup ,p) - - w l i^Ttrm 

v t/i./t/ 2 ,/fe(k + q))-/fe(q)) 
+ M &(k + q)-&(q) ' ^ ' 

/ is the Fermi distribution of £,j (q) = e q — fij, k' = p + p' and k = p — p', where p, 
are the momenta of the two pairing spin states. Induced interactions due to particle- 
hole loop diagrams from j = 3 , V spin states are responsible for the second sum 
in Eq. d44l and has opposite sign. For two components, however, it is absent. 

In the two-component spin-balanced system the induced interactions effectively 
leads to a second order correction to the scattering length a — > a + (2/3n)\n(4e)kpa 2 
ll50ll5D . where a = ayi and Un = Uq = 47th 2 a/m. Including this correction in Eq. 
d43l reduces the gap by a factor (4e) 1,/3 ~ 2.2 



fix 1 " 


n 


[ - 1 E F exp 


2akf 



(45) 

The induced interactions consists of a repulsive direct part and an attractive part 
due a loop diagram which is therefore proportional to the number of components v. 
The induced interactions therefore scale as (3 — v) in a system of V spin balanced 
multi-components with the same scattering length. The corresponding gap is 

A = (4e) v/3 - l A . (46) 

For spin polarized systems or for the 6 Li systems with three different scattering 
lengths the pairing gap is described in Refs. l52l . Adding bosons enhance pairing 

Ell- 
in the unitary limit kp \a\ xLl the gap is of order the Fermi energy IflOl . Extrap- 
olating (|45l l to ak F — > ±°o J3] gives a number A = 0.49£> close to that found 
from odd-even staggering binding energies A = 0.54Ep calculated by Monte Carlo 
fl4l . The crossover model of Leggett described below gives a somewhat larger gap 
A Leggett /Ep = 0.69. Such values of order the Fermi energy are one or two orders 
of magnitude larger those found in metals and high temperature superconductivity 
and was met by disbelief among condensed matter physicist and the crossover was 
considered academic. Only after the realization and confirmation of the BCS-BEC 
crossover in experiments with ultracold atoms was its significance acknowledged. 



2 Inquiring into details about their calculation some 40 years later Gorkov only remembered "...that 
it was a particular difficult calculation!" 
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Fig. 5 Superfluid transition temperatures from the Leggett crossover, Gorkov and NSR. 



It also follows from the BCS gap equation that the gap and thus superfluidity and 
superconductivity all vanish at a critical temperature 



T c = —A -0.5674. 

K 



(47) 



for any weak interaction (Ce = 0.577... is Euler's constant). In the unitarity and 
BEC limits the critical temperature is no longer proportional to the gap. 



3.3 Leggett's BCS-BEC Crossover 

Eagles and Leggett made the important step of connecting BCS pairing to a BEC 
of molecular bosons via crossover. Eagles gave an early treatment of the crossover 
from BCS superconductivity to a BEC in the context of systems with low carrier 
concentrations (9]. For pedagogical reasons we will follow Leggett who solved the 
gap equation of Eq. (|42T > and number conservation 



N - 



E 

k 



1 



£k-M 



(48) 
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self consistently in the crossover. The beauty of Leggett's crossover model Eqs. ( |42] | 
and d48l > is that it describes the crossover continuously as function of the parame- 
ter x = l/akf from the BCS limit (x —> — °°) over the unitarity limit (x = 0) to the 
BEC limit (x — ^ +°°). It is important to include (pairing) interactions in the chemical 
potential through number conservation as it changes the chemical potential dramat- 
ically. 

In the dilute (BCS) limit the gap equation leads to the standard BCS gap of Eq. 
d43l ) - not including the Gorkov correction. The chemical potential is jj. = E F and 
does not include the standard mean field Hartree-Fock correction of a dilute gas. 

In the BEC limit the pairing gap approaches 



A =AE F /yj3%ak F . (49) 

The chemical potential approaches half of the molecular binding energy Em = 

—h 2 /ma 2 , 

h 2 7th 2 na 

M = -^^ + , (50) 

plus the BEC mean field corresponding to a molecular scattering length of om = 2a. 
Four-body [54], Monte Carlo calculations |[T5ll and experiments lfl6l do, however, 
indicate that «m — 0.6a 

On the BCS side the minimum quasiparticle energy is A and occur when k = k F . 
On the BEC side the chemical potential is negative and the minimum quasiparticle 
excitation energy is the quasiparticle energy for k = 



E qp = + (51) 

The quasiparticle energy is observed in the spin excitation response function IT761 
157). 

In the Leggett model the gap vanishes in the BEC limit at a critical temperature 
jdissoc _ |e M | j /l n (|e M |/J5 Jf ) 3 / 2 ; which is smaller than the quasi-particle excitation 
energy. It is not the on-set temperature for superfluidity but rather a molecular pair 
dissociation temperature 0111 . The onset of superfluidity occurs at a lower critical 
temperature for a molecular BEC, Tf EC . 

The Leggett model fails in the BEC limit because of the basic assumption that 
only opposite momenta fermions (zero total momentum pairs) can pair. Pairs with 
non-zero momenta are thermal excitations of molecular bosons and including such 
degrees of freedom lowers the critical temperatures increasingly towards the (molec- 
ular) BEC limit. 
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3.4 NSR 



The model of Nozieres & Schmit-Rink (NSR) ifTTI extends the Leggett model so 
that it correctly describes the critical temperature in the molecular BEC limit. We 
will therefore give a brief outline of NSR with emphasis on how pair motion is 
included and how it corrects the Leggett model in the BEC limit. Also, the NSR 
approach is formulated such that it applies to optical lattices discussed in sec. (0 
with few but crucial differences. 

In the NSR model the (molecular) pair momentum q is included in the two- 
particle correlation function. To lowest order it is given by the propagator for two 
free atoms 

T 

n(q,co v ) = 77 E Go(k,ico m )G (q-Ki®v - ico,„) 
k,a>,„ 

1 -/( £ q/2+k) -/(£q/2-k) ^ 
V k COv + 2/X- £ q /2+k - £q/2-k ' 

In Eq. (TS2b the Matsubara frequencies o,„ = mlnTi have been summed over integers 
m. 

We now scale the potential by A and sum up interactions energies (UTI) n from 
n = 1,2, ... ladders. We obtain for the expectation value of the interaction energy in 
state X gives 

From the Hellmann-Feynmann theorem we now obtain the interaction part of the 
free energy 

flj* = / i- (Atf) = T £ ln[l - C/n(q, my)] . (54) 

As in the gap equation d42b large momenta contributions are removed by replacing 
UIl by UoTIr, where the renormalized propagator is 

n r (q,© v )=n(q,«v) + ^E^r. (55) 

Inserting into Eq. d54"t gives the thermodynamic potential 

i2=i2o + r£ln[l-f/ n r (q,fi) v )], (56) 

q.fflv 

where Qq = — 2T^ q ln[l +e(% _ M)/ 7 '] j s the free energy for non-interacting Fermions. 
The frequency sum in Eq. 061 ) can, using the residue theorem, be converted to an 
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©-integral of the Bose distribution function (exp(o/r) — 1) around the real axis, 
where the logarithm in Eq. d56l ) has a cut. 

In the BEC limit, where n is large and negative, it follows from Eq. d52l that 
1 — Uon r (q, (o) is proportional to co — £m(#) + 2/1, where 



h 2 n 2 q 2 

M? = J + -p-, (57) 

ma L 4m 



is the molecular binding and kinetic energy of a pair with mass 2m and momentum 
q. Therefore the derivative of ln[l — U(,TI r (({, (0 V )} with respect to jj, has a pole at 
co = £m(<?) — 2/j. with residue 2. Using the residue theorem again we obtain the 
number equation 

N ~~lijl - 2 ^ e (e M M-2^)/T Z i ' (58) 

Here we have ignored the contribution to the thermodynamic potential Qq from free 
Fermions which is negligible in the BEC limit. The number equation is simply that 
of a free Bose gas as opposed to that of a Fermi gas Eq. d48l ) in the BCS limit. 

Whereas the number equation changes qualitatively from free Fermi atoms to 
free Bose molecules in the crossover, the gap equation is unchanged and given by 
Eq. (142) . In the BEC limit the gap equation simply yields that the chemical potential 
is half the molecular binding energy, jU = —Ti 2 /2ma 2 . The number equation now 
gives 

K h 2 n 2 / 3 

— ~0.218£ F , (59) 



[2f(3/2)] 2 /3 m 

i.e., the critical temperature is independent of the pairing interaction and given by 
the T C BEC for the molecular BEC. At T C BEC all bosons are thermally excited with 
none remaining at zero momentum. 

In general the gap equation ( 1381 and number conservation N = d£2 /dfl with 
£2 = £2q + Cli nt from Eq. (T54b have to be solved self-consistently for the critical 
temperature and the chemical potential in the BCS-BEC crossover. The result for 
T c is shown in Fig. (0. One notices that NSR predicts a maximum near the unitar- 
ity limit. Recent experiments find T c /Ep = 0.19 ±0.2 l20l whereas earlier found 
T c /Ef = 0.29 ± 0.3 J4). This is compatible with the gap A = 0.44£> found in 1531 . 

Although the NSR model is a qualitative improvement of the Leggett model by 
correctly describing T c in the BEC limit, it still omits a number of effects such as 
particle-hole contributions, selfenergies and induced interactions from the medium. 
Consequently, the Hartree field and Gorkov corrections are not included, and in the 
BEC limit the molecule-molecule scattering length 2a implied by Eq. d50l is a factor 
~ 3 too large 1541 . More elaborate models 156113711581 include some of these effects 
and do not find a maximum for T c /Ep around the unitarity limit. 
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4 Atomic traps and Nuclei 

Traps are necessary for confining the atoms and creates a density distribution that is 
maximal at central density and decrease towards the cloud size R where it vanishes. 
Thus all densities are present at once which can be difficult to separate experimen- 
tally in order to extract the (n,T) phase diagram. 

From a nuclear point of view harmonic oscillator traps are wonderful toy systems 
since they provide confined systems, which are a first approximation to the nuclear 
mean field potential. Additionally we can tune the interactions and therefore the 
level splitting and pairing. Also millions of Fermi atoms can be confined and not 
just the ^,250 neutron and protons in nuclei limited by fission. Therefore one can 
study the crossover from few to infinite number of particles which is necessary in 
order to e.g. link pairing in nuclei to that in neutron and nuclear matter. Similarities 
between the (prolate) nuclear mean field and cigar shaped optical traps can also be 
exploited in studies of shell structures, pairing and collective modes. 
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4.1 Scaling in atomic traps 

In a noninteracting two spin balanced Fermi gas in a spherically symmetric HO 
potential at zero temperature the particle fill HO shells of energy nTico, n= 1,2, ..,Hf 
up to the Fermi shell np = (3./V) 1 / 3 with Fermi energy Ef = iifTiCO. For a sufficiently 
large number N of particles confined in a (shallow) trap the system size R is so long 
that density variations and the extent of possible phase transition interfaces can be 
ignored and one can apply the local density and the Thomas-Fermi approximation. 
Here the total chemical potential is given by the sum of the harmonic oscillator (HO) 
trap potential and the local chemical potential /i, = dE / dNj 



which must be constant over the lattice for all components i = 1 , 2 , . . . It can therefore 
be set to its value at the its edge /?,-, which gives the r.h.s. in Eq. d60b . We shall 
mainly discuss spin balanced two component systems where the chemical potential 
and radii are equal. The equation of state determines jll(n) in terms of the universal 
function of Eq. ( l60l >. 

In the dilute limit and in the unitarity limits ji = £,Trkp/2m, where § = 1 
and | = 1 +/3(0) respectively. In both cases Eq. (f60b gives n(r) = kp{r)/3n 2 = 
n o(l — r 2 /R 2 ) 3 / 2 , where no = ^ _3//4 (2«/?) 3 / 2 /3n 2 a\ sc is the central density and 
R = t, X I A '\j2nFa osc the cloud size; a osc = ^fi/mco is the oscillator length. The at- 
traction contracts the gas to a higher central density such that R °< kp (0) t, X I A '. 
The total energy of the trapped gas is E/N = {3/%)^l 2 n F hco. 

For a more general equation of state P °c or fj « «'' with poly tropic index y 
the equilibrium density is n eq = «o(l — r 2 /R 2 ) l /Y, where R 2 — 2(y+ l)Po/ynoma> 2 . 
Po and «o are the pressure and density in the centre of the trap. In both the non- 
interacting and the unitary limit y = 2/3 whereas y = 1 in the BEC limit and for a 
bose gas. 

As suggested in lfl2l it is convenient to measure the density along the axial axis 
integrated over transverse cross section n(z) = f n{r)dxdy. From the Gibbs-Duhem 
equation (dP/dji)T = n using Eq. d60t we obtain the pressure in the centre along 
the axial axis by integrating over transverse coordinates 



At the same time the chemical potential ji(z) =m(O z (R 2 — z 2 ) is known by measuring 
the size R(z) of the cloud along the z-axis. Thus P(jtt) can be measured and the 
equation of state extracted at any temperature. 

In a Fermi liquid we obtain from Eqs. (0 and (fl4] i that the pressure w.r.t. that in 
a free Fermi gas is 




m(0 2 R 2 



i ■ 



(60) 



57T 2 m* T 2 

P( X ,T)/P FG = 1 + P(x) + - -J. 

o m E F 



(62) 



Crossovers in Unitary Fermi Systems 



27 



Detailed measurements of the intensive variables (T,P,ji) gives J3(0) = —0.49(2) 
and m*/m= 1.13(3) El. 



4.2 Collective modes 

Tickling the trapped atoms sets them into oscillations at certain eigen-frequencies 
called collective modes. Such giant dipole and quadrupole modes have been im- 
portant for studying nuclei |59l . The collective modes can be calculated from the 
equation of state and the Euler equation 

mn-^- = -V(P + nmY^a)frf). (63) 

i 

where v = dr/dt is the local velocity. The last term is the gradient of a generally 
deformed HO potential. The Euler equation can be solved analytically for polytropic 
equation of states. For spherical symmetric 3D traps one finds collective modes at 
eigen frequencies aty for a mode with rj = 0, 1,2, ... nodes and angular momentum 
/ given by l60l 

toll 

-JE=l+2 n y tj+/ + 1/2) + 1 . (64) 

Similar expressions exist for modes in deformed traps ll6D which have been mea- 
sured in detail BlfTTl. 

For attractive interactions the system is superfiuid at zero temperature and there- 
fore irrotational with a quenched moment of inertia. As the temperature is increased 
above T c the moment of inertia increase to rigid value and hereby the critical tem- 
perature T c = 0.19(1)£V is found 1621 . Deformed system can be rotated but has a 
quenched moment of inertia / = 5 2 /, (?(£ / l63l . 



4.3 Shell structure and pairing in atomic traps 

The pairing gap generally increase with the level degeneracy or density of states 
|59l . For pairing in a single level of angular momentum / or shell rip the pairing gap 
scales with the number of degenerate states (21 + 1) or (rip + l)(«f +2) respectively. 
If levels are split or distributed the gap depends on the density of levels around the 
Fermi level as in the uniform case of Eq. (143) . 

For very weak interactions the level splitting of single particle levels l=riF,np — 
2, 1 or in the HO shell uf is smaller than the pairing gap when there are suffi- 
ciently few («/r^S10) atoms in a HO potential. The pairing occurs between all states 
in the shell leading to the supergap l65l 
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Fig. 7 Diagram displaying the regimes for the various pairing mechanisms (see text) at zero tem- 
perature in h.o. traps vs. the number of particles N = rtp/3 and the interaction strength a. The 
dotted lines indicate the transitions between single-shell pairing A = G, multi-level, single-level, 
and multi-shell pairing. At the dashed line determined by 2Gln(yjz/r) = ha> the pairing gap is 
A ~ hco, and it marks the transition from multi-shell pairing to bulk superfluidity Eq. (45). The 
pairing gap is A = 0.54-Ef above the full line p|a| 3 > 1, which separates the dilute from the dense 
gas. From 1651 



G= 3V2^+3 j*^ (65) 
\5k z a osc 

For more particles or stronger interactions the situation becomes more complicated 
depending on the gap size with respect to the level splitting as shown in Fig. (0. 
Increasing the number of particles cause level splitting so that pairing is reduced to 
multi-level and eventually single level pairing. Increasing the interaction strength 
increase pairing to nearby shells ll65l referred to as multi-shell pairing in Fig. (|7]). 
The general case can be solved within the Bogoliubov-deGennes equations ll65l[3TI . 
For stronger interactions and many particles the pairing approaches that in a uniform 
or bulk system, and eventually dense system with \a\k,F > 1 and the unitarity limit 
is reached. 

The various pairing mechanisms and phases of Fig. (0 can be studied at low 
temperatures is which hopefully will be reached in the near future. 
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4.4 Pairing in nuclei 

We can exploit some interesting similarities between pairing in nuclei 1591 64) and 
that of Fermi atoms in traps l65l |3T1 . To a first approximation the nuclear mean 
field is often taken as harmonic oscillator (HO) potential just as the optical traps. 
Secondly, the residual pairing interaction between nucleons is taken as a short range 
(delta function) interaction as for cold atoms. 

Before these result can be applied to nuclei there are, however, a number of dif- 
ferences that must be taken into account. Large nuclei have approximately constant 
central density po ~ 0.14 fm~ 3 and Fermi energy Ep in bulk. Therefore the HO 
frequency, which is fitted to the nuclear mean field, decreases with the number of 
nucleons A = N + Z, where N now is the number of neutrons and Z the number of 
protons in the nucleus, as Tico ~ Ep/np ~ 41M6V/A 1 / 3 . In the valley of ^3 -stability 
the number of protons is Z ~ A/(2 + 0.0155A 2 / 3 ). Therefore in medium and large 
nuclei the difference between the Fermi energies of protons and neutrons exceeds 
the pairing gap so that pairing between protons and neutrons does not occur. 

Secondly, the nuclear mean field deviates from a HO potential by being almost 
constant inside the nucleus and vanish outside. The resulting net anharmonic nuclear 
field is stronger and opposite in sign to the corresponding (anharmonic) mean field 
in atomic traps. Therefore, the level splitting is larger and the ordering of the 1- 
levels is reversed. In addition, a strong spin-orbit force splits the single particle 
states of total angular momentum j = Z± 1/2, such that the j = np + 1/2 is lowered 
down to the shell (np — l) below. The level splitting can be parametrized by a single 
parameter taken from analyses of nuclear spectra. It increases with shell number 
up to «f ~ 6 for heavy nuclei. Due to the strong spin-orbit force the j = I ± 1/2 
states are split and the j = np + 1/2 is lowered down to the shell below. The magic 
numbers become N,Z = 8, 14,28,50,82, 126, 184, .., etc. rather than the h.o. filled 
shell particle numbers N,Z = 2,8,20,40,70, 1 12, 168,240, etc. 

The pairing gaps and quasi-particle energies can now be calculated by solving the 
Bogoliubov-deGennes gap equation l65l which are shown in Fig. (0 for neutrons. 
The strong level splitting in nuclei has the effect that pairing is strongest when the 
shell is half filled and weak near closed shells simply because there are fewer states 
available for pairing, i.e., the level density is smaller. Averaging over several shells, 
however, the mean gap is well approximated by the supergap because the reduction 
in pairing due to level splitting is compensated by additional pairing to nearby shells. 
Since Tico scale as ~ A -1 / 3 and a osc <x n l J 2 the single-shell pairing gap also scales 
as G ~ A -1 / 3 . Therefore, the pairing gaps in light and medium mass nuclei scale 
approximately as OT1 



A ~G~ 



\a\ 5.5MeV 
0.41fm A'/3 



(66) 



As shown in Fig. ([8]) the supergap does not depend on the level-splitting and is 
therefore a robust prediction for the average magnitude and mass scaling of pairing 
gaps in nuclei. 
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Fig. 8 Neutron pairing energies vs. the number of neutrons. The experimental odd-even staggering 
energies A' 3 '(N) are averaged over isotopes |67] [66| . The calculated gaps A and quasi-particle 
energies E are obtained from the gap equation (see text) with effective coupling strength a = 
—0.41 fm. The supergap G is shown with dashed line. See [311 for details and corresponding plot 
for protons. 

The data on neutron and proton pairing is obtained from the odd-even staggering 
of nuclear binding energies B(N,Z). It has been shown that mean field contribu- 
tions can be removed ||66l by using the three-point filter A^(N) = (—l) N [B(N — 
l,Z)/2 + B(N+l,Z)/2-B(N,Z)]. We compare in Fig. © to the experimental 
A ( 3 ) (N) averaged over isotopes with the calculated gaps. The analogous for protons 
A^(Z) averaged over isotones compare similarly to Bogolibov-deGennes calcula- 
tions 16511311 . In the calculations the effective coupling is the only adjustable pa- 
rameter which is fitted to the experimental data. For both neutrons and protons we 
extract a ~ —0.41 fm. 

Considering the simplicity of the model it describes a large number of experi- 
mental gaps fairly well on average. In a number of cases, however, the calculated 
pairing gaps differ significantly from the measured neutron gaps. Some of these 
deviations can be attributed to the crude single particle level spectra assumed. 

The pairing in nuclear matter can now be estimated once the effective interaction 
has been determined using Eq. d43l ). Note that the effective interaction includes in- 
duced interactions. Inserting a = —0.41 fm and kf = 1.3 fm -1 at nuclear saturation 
density, po = 0. 15 fm~ 3 , we obtain the proton and neutron pairing gaps 

4~l.lMeV, (67) 

in the bulk of very large nuclei and in symmetric nuclear matter at nuclear saturation 
density. This number is compatible with earlier calculations l39l of the 'So pairing 
gap in nuclear and neutron star matter around normal nuclear matter densities. 
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Neutron star matter has a wide range of densities and is asymmetric, Z/A ^0.1, 
above normal nuclear matter densities. One can attempt to estimate of the pairing 
gaps as function of density from the gap in bulk, Eq. d43l l, with a ~ —0.41 fm 
and the neutron or proton Fermi wave numbers kp' Z = {3n 2 Pn ,z) 1 ^ as function 
of densities. However, the effective interaction a is density dependent. At higher 
densities we expect that the effective interaction becomes repulsive as is the case for 
the nuclear mean field at a few times nuclear saturation density. At lower densities 
the effective scattering length should approach that in vacuum which for neutron- 
neutron scattering is a^So) ~ — 18 fm. 

In spite of several simplifying approximations in this atomic trap model for nu- 
clear pairing it provides at least a qualitative description of pairing in most nu- 
clei. The effective interaction a = — QAlfm is the same for both neutrons and pro- 
tons which reflects that both are close to forming bound states. At central densities 
x = 1 /akf ~ —1.9, i.e. on the BCS side of the unitarity limit. 



4.5 Quark and gluon matter 

The Coulomb and QCD interactions ~ g 2 /r or their Fourier transform ~ g 2 /q 2 are 
long range and therefore contrary to the short range interactions underlying univer- 
sality. The long range QED and QCD interactions requires screening of infrared 
divergences ||68l whereas in unitary gases ultraviolet cutoffs are provided by renor- 
malizing the short range interaction in terms of the scattering length. Never the less 
a number of similarities have appeared. 

Confinement is caused by strong color fields between quarks and gluons which 
forms bound state hadrons at temperatures below T^JTqgp — 160MeV and densi- 
ties below a few times normal nuclear matter density, po = 0.15fm _1 . These may 
again form a nuclear liquid at temperatures below the nucleon gas critical tempera- 
ture T<Tnq ~ 15MeV for a four component nuclear system as discussed above. In 
the hadron gas the interparticle distance is larger than the confinement interaction 
range, and the quark molecules can therefore be viewed as a molecular gas. How- 
ever, color neutrality requires either Fermi molecules of three quarks (baryons) or 
quark-antiquark Bose pairs (mesons). Of the two phase transitions in the molecu- 
lar gas of hadrons only the neutron gas can be viewed as a BCS-BEC crossover as 
discussed in sec. (12.6b . The instability and first order transition of a four-component 
nuclear gas was discussed in sec. ( 12.61 l can be viewed as multi-component crossover. 
Here the effective scattering length is that between neutrons and protons a ~ 20fm, 
i.e. on the molecular BEC side, and the nuclear matter is in the unitarity limit in 
the sense that x = 1/ akf ~ + although they now are kept in place by short range 
repulsive forces. The transition from nuclear to quark matter or a hadron gas to 
a quark-gluon plasma is quite different because the interaction range is always of 
order the interparticle distance. The strong short range repulsive forces between nu- 
cleons leads to a strongly correlated wavefunction as in hard sphere scattering and 
the strongly correlated nuclear fluid. In this case the Gutzwiller approximation of 
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Ref. (46j is valid and the matter does not undergo a ferromagnetic transition sepa- 
rating the unlike components (see Fig. ((TJ). As seen from the quark-gluon side one 
can define a "crossover parameter" x ~ Aqcd /kf which at high densities is in the 
unitarity limit. With decreasing densities the running coupling constant Aqcd di- 
verges when the quark-gluon plasma undergoes a first order transition to a hadron 
gas with x^jl. 

Quark pairing also have similarities to pairing in nuclear and multi-component 
atomic gases. Color superconductivity between the 2 spin, 8 color and 2-3 flavor 
states is very sensitive to flavor imbalance JS) . If the strange quark mass is small the 
up, down and strange quark Fermi levels are close. If pairing is sufficiently strong 
such that the gap exceeds the Fermi level splitting all flavors can pair. It is amazing 
that properties of pairing in quark matter, which may exist in unaccessible cores of 
neutron stars, can be studied in tabletop experiments with ultracold multicomponent 
spin-imbalanced Fermi atomic gases. 

Another interesting similarity is elliptic flow. The overlap zone in semi-central 
high energy nuclear collisions is prolate (cigar shaped) as are optical traps. In sub- 
sequent expansion the hydrodynamics forces stronger expansion in the direction 
where systems is narrowest initially. This makes the momentum distribution az- 
imuthally asymmetric - referred to as elliptic flow. Results indicate that both "liq- 
uids", the QGP/hadronic and ultracold unitary atomic gas, expand as almost perfect 
fluids initially i.e. their viscosities are record breaking low |69l . 

Correlations can reveal the quantum phase structure. Originally Hanbury-Brown 
& Twiss measured the Bose-Einstein correlations between stellar photons and de- 
termined the diameters of nearby stars. Similar correlations between mesons in high 
energy nuclear collisions have been exploited to determine the freezeout size of 
the collision zones |70|. Also Fermi anticorrelations have been observed between 
baryons. In ultracold atomic systems analogous "noise" correlations have been 
found near the BCS-BEC crossover due to pairing JT5J. Bragg peaks have been 
observed for bosons in 3D l7T| and 2D 1721 lattices, and dips for 3D fermions in 
El. 



5 Optical Lattices 

A standing wave can be generated by shining lasers onto traps which generates a 
3D cubic or 2D square lattice potential with spacing A ~2n/k, which is half the 
laser wavelength. The lattice height Vq varies with laser intensity and the scattering 
length. In the tight binding approximation the onsite coupling 




(68) 



and hopping parameter 



t = E R (2/^ 3 e- 2 ^, 



(69) 




where Er = h 2 k 2 /2?n is the recoil energy and £, = (Vq/Er) 1 / 4 . Varying Vo and the 
scattering length a near Feshbach resonances allows us to tune the Hubbard param- 
eters U and t . 

The lattice constant A introduces another lengths scale so that the thermodynamic 
quantities generally depend on both density and interaction strength. Therefore uni- 
versality is broken except for low filling which corresponds to dilute Fermi gases. 

The Hubbard Hamiltonian on a D-dimensional lattice describes optical lattices 
that are sufficiently deep for a one-band model to apply, 

H = U o,a'nion ia i -t £ a] a a ja . (70) 

!,(J<(j' 

Here a ia is the Fermi creation operator of the hyperspin states a = 1 , 2, . . . , V, n iCT = 
a^fl/CT the density and (ij) denotes nearest neighbours with hopping parameter t . 
Due to particle-hole symmetry results for a given site filling n also applies to 2 — n. 

The Hubbard model is of fundamental importance in condensed matter physics 
where it explains properties of Mott insulators, antiferromagnets, d-wave supercon- 
ductivity, etc. The Hubbard model can only be solved in ID for fermions but a num- 
ber of models and numerical calculations have given insight in the various phases 
and order parameters. Optical lattices may tell us the answers in a few years. 



5.1 On-site attraction and Pairing 

We start by investigating s-wave pairing between to spin balanced states due to an 
attractive on-site interaction U < 0. The mean field gap equation for singlet super- 
fluidity at zero temperature is l74l 
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Fig. 10 Superfluid transition temperatures in a 3D lattice at half filling for the NSR model, the 
Heisenberg model within (Weiss) mean field T c = 1.57 and Monte Carlo. "Gorkov" is the weakly 
attractive limit of Eq. {74) reduced by a factor 3.9 due to induced interactions. 



U ,, s tanh(£„/2r) 

T XU(P\P)\ L P/ \ (71) 



where M is the number of lattice points and E p = yj (e p — fj,) 2 + A 2 with e p = 

2fEi=i^)(l-cospi)- 

The density n = 1 — ^ p (ep/£ , p)tanh(£ , p /2r)/M, is also the filling fraction in 
units where the lattice constant is unity. Momenta are in the first Brillouin zone only 
\pi\ < % taking A = 1 for convenience. 

Due to particle-hole symmetry results also apply replacing the density by (2 — n) 
and chemical potentials by (4Dt — fl), where 4Dt is the bandwidth in D dimensions. 
At low filling e p ~tp 2 and results also apply to uniform systems replacing t = 1 /2m. 



5.1.1 Pairing in the dilute limit 

We first calculate the zero temperature gaps for weak attraction \U\ <C t excluding 
induced interactions. The s-wave gap A p > = Aq is then momentum-independent be- 
cause the interaction U is. The gap equation (|7TTi reduces to: 1 = —(Uo/M)Y, v 1 /2£p, 
from which the gap can be calculated analytically at low and half filling to leading 
orders in 3D and 2D (for ID see, e.g. l75l ). 
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The 3D lattice has a critical coupling ||76l 



U e = -M/(£l/2e q ) = -8v/2r/[£/^ 

q 1=0 



l)!!/2 2/ /!] ~ -7.913f,(72) 



where a two-body bound state can be formed at zero density (jj. = 0). This threshold 
corresponds to the unitarity limit of infinite scattering length in the uniform system 
and naturally enters the sum in the gap equation. The dilute gap becomes to leading 
orders 



[3D 



r^exp 



'4n 2 t 






[ k F 




-k)] 



,n<l, 



(73) 



at low densities where /I = tkp <?C t . The level density at the Fermi surface N(jJ.) = 
kf /{An 2 t) enters in the exponent as in standard BCS theory. U c acts as a cutoff in 
the one-band Hubbard model which in the uniform (continuum) limit is absorbed 
into the scattering length a. Thus Eq. d73l is the finite lattice equivalent of the gap in 



the uniform system Aq D = (8/e 2 )iiexp(7t/2k F a), with U~ l — U^ 1 = t/ ~\ where 
Uq = Ana/m as in the uniform case, such that threshold U = U c corresponds to 
\a\ = oo. 

Near half filling the level density is surprisingly constant N(n) — 0.143/f in a 
wide range At < jj, < 8f around half filling. To leading orders the gap becomes ||74l 



A 3D 



at exp 



1 



1. 



(74) 



where the prefactor a ~ 6.544 can be calculated numerically. 

The 2D lattice has a superfiuid s-wave gap that be calculated from the gap equa- 
tion El 



[2D 



y/STtjlt exp 



Am 



u 



,n< 1, 



(75) 



when the density n = k 2 F /2n is small. The 2D level density is N(n ) = 1/ {Ant). Eq. 
( TTSb assumes that the on-site coupling and the gap are small such that Aq D ^[1. This 
is not fulfilled at sufficiently low densities where instead Aq D = Ant exp(8nt /U), 
which also is the two-body binding energy in 2D at zero density. There is always 
a two-body bound state in 2D with purely attractive interaction and therefore U c 
vanishes. In 3D a similar pair condensate (a molecular BEC) appears when a > 
corresponding to U < U c . 

For intermediate fillings the pairing is more complicated in 2D. Near half filling 
n~ 1 the level density N(s) = ln(16f/|e — At\)/(2n 2 t) has a logarithmic singularity 
due to the van Hove singularity. Calculating the r.h.s. of the gap equation there- 
fore gives a double log: 1 = \U\/(Ak 2 ) In 2 (32t/AQ D ), to leading logarithmic orders, 
resulting in the gap when n ~ 1 lf77l 
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(76) 



5.1.2 Induced interactions on the lattice 

In 3D the induced interactions are at low densities n< 1 the same as in the uniform 
case Eq. ( PPfl i and reduce gaps and T c by a factor (4e) 1 ^ ~ 2.2. Near half filling 
the reduction factor is ^3.9 for weak interactions. For stronger interactions the 
induced interactions are enhanced J52j leading to a reduction of the gap by a factor 
~ 25 when U = -3t. 

In 2D the induced interactions suppress the gap by a factor e in the dilute limit 
l78l but are divergent at half filling due to the singular level density ll52l . The in- 
duced interactions have been calculated for general spin polarization and number of 
components in 3D and 2D optical lattices (53). 



5.1.3 Crossover and the BEC limit 

As the interactions become stronger the Fermi gas gradually undergo crossover to 
a molecular BEC. On the lattice, however, the finite bandwidth automatically pro- 
vides a cutoff which affects the crossover and changes the BEC limit on the lattice 
qualitatively from the uniform case. We shall describe this crossover employing the 
NSR model as above for the uniform case but with a finite bandwidth. Our starting 
point is again the two-particle correlation function of Eq. d52l . In the BEC limit 
(large negative chemical potential) we find from Eq. d52l that 1 —UIl(q, Co) is pro- 
portional to co — £m(i) +2/X as in NSR, where now the molecular boson energy 
is 

e M (q) = -U+jY, [1 -cos(g,)] (77) 

'=1,3 

with / = 4t 2 /U. Therefore the number equation reduces to the condition for BEC 
in an ideal Bose gas as Eq. d58l l 

N ~ ^7 J -~ 2 E e ( £M (q)-2 M )/r_ 1 ' (?8) 

but with the molecular boson energy of Eq. ( 1771 ) in stead of d57l ). The gap equation 
at T c yields a chemical potential fi = —U /2 in the BEC limit, and inserting this in 
the number equation gives a critical temperature 



(79) 



Crossovers in Unitary Fermi Systems 



37 



when n<l. For n = 1 we obtain T c = 1.717. 

The crossover in optical lattices is shown in Fig. ( [Tol l at half filling. The critical 
temperature decreases with interaction strength in the BEC limit and thus differs 
from uniform systems where it approaches a constant Tf EC . 

The Hubbard model can in the BEC limit be mapped onto the Heisenberg spin 
model in a magnetic field for which accurate calculations exist, which allows us 
to check the NSR and other models ll79l . Half filling corresponds precisely to zero 
magnetic field, i.e., 

H = j£S r Sj, (80) 

(y> 

and the Neel temperature corresponds to Tn = T c in the BEC limit. In (Weiss molec- 
ular) mean field theory Tn = 1 .57, quantum Monte Carlo Tn = 0.957, high tempera- 
ture expansions Tn = 0.907, as compared to T c = 1.717 in NSR. Extending the NSR 
model with charge and spin density fluctuations gives a smaller value T c ~ 0.47 
because molecular repulsion is overestimated |79i|. 



5.2 On-site repulsion 

On-site repulsive interactions U > generally disfavours doubly occupied sites and 
lead to a Mott gap at half filling of order ~ JJ jl at very low temperatures. Mott 
insulator (MI) transitions are observed for bosons in traps at fillings n = 1 , 2, 3,4,5, .. 
lfl6l . For fermions the phase diagram is more complex, an antiferromagnet (AF) 
at and near half filling, a paramagnet (PM) otherwise and possible a ferromagnet 
for very strong repulsion l86l . Furthermore the AF phase may be unstable towards 
phase separation to two coexisting phases: an AF at half filling and a PM at filling 
slightly below or above half filling. The co-existing phase may also be unstable 
towards stripes. In 2D these phase also compete with d-wave superfluidity. 

5.2.1 Antiferromagnetism 



At half filling fermions are known to form a MI in ID [80] whereas in two and 
higher dimensions an antiferromagnetic (AF) insulator is found in Monte-Carlo cal- 
culations. The AF alternating spin order on the lattice is driven by the simple fact 
that hopping can occur to a neighbouring site only if it is occupied by an opposite 
spin, which generates a super-exchange coupling. 

We will here study the AF phase and its transition to a paramagnet at the Neel 
temperature 7]y. Remarkably, it undergoes a crossover very similar to the critical 
temperature described above for attractive interactions, i.e. Tn ~ T c . For weak re- 
pulsion Tn can be determined from the mean field gap equation 
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U y tanh(£ p /27]y) 



(81) 



which is identical to that for T c of Eq. (ITTl i except for the sign of f/. Therefore Tn = 
T c within mean field. Induced interactions can be included as for T c as described 
above and therefore reduce Tn by a factor ^3.7 at half filling IBTl which is very 
similar to the reduction of T c by induced interactions as discussed above. 

For stronger coupling fluctuations reduce the mean field gap and Tn as for the 
NSR model of T c . At very strong coupling the repulsive Hubbard model can at half 
filling also be mapped onto the Heisenberg model with coupling J. Therefore Tn 
is also in this limit given by T c as function of \U\, and the details of the crossover 
ll82l[83l are quite similar to the BCS-BEC crossover. 

The AF phase transition masks a possible MI phase. By frustrating the system 
one can, however, remove the AF order and observe the MI transition by a vanishing 
conductivity or compressibility. Another way to suppress the AF order is to include 
multi-components which naturally upsets the alternating spin up/down order in a 
two-component AF. The limit of many components actually correspond to a Bose 
system where MI transitions are found at every integer filling. 

5.2.2 Phase separation 

Mean field theory provides a first impression of the phases competing for the ground 
state and has the advantage that it is computationally simple as compared to more 
complicated theories. The MF equations for the Hubbard model are standard and 
we refer to, e.g., Refs. 07411841 . The energy densities can be calculated within the 
Hartree-Fock approximation for the paramagnetic (PM), ferromagnetic (FM), anti- 
ferromagnetic (AF) and phase separation (PS). 

At low density n<l the ground state is that of a dilute paramagnetic (PM) gas 
with energy 



In units where the lattice spacing is unity (A = 1) this energy per site is also the 
energy density, and the density is the site filling fraction. 

Near half filling Epm = —{4/n) 2 t + U/4, which becomes positive when the re- 
pulsive interaction exceeds U/t > 64- /n 2 ~ 6.5. The PM phase is then no longer the 
ground state. In a state with only one spin the antisymmetry of the wavefunction 
automatically removes double occupancy and the repulsive term Un 2 /4 in Eq. d82l 
disappears. Such a ferromagnetic (FM) state always has negative energy Efm < for 
n < 1 and is a candidate for the ground state. AF, linear AF [83] and stripe phases are 
other competing candidates. Furthermore phase separation (PS) and mixed phases 
can occur near half filling. When U^Jt the ground state of the 2D MF Hubbard 
model undergoes transitions from an AF at half filling to a mixed AF+PM phase 




(82) 
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for doping 8 = 1 — n up to a finite (U dependent) value |5| < 8 S where after a pure 
PM phase takes over. |84l For larger U the phase diagram is more complicated with 
a pure as well as mixed FM phases between the AF and PM phases. A finite next 
neighbour hopping term t ' makes the phase diagram asymmetric around half filling, 
extends the AF phase and changes the phase diagram considerably. 

At half filling n = 1 the ground state is an AF. Near half filling < |<5| <C 1 the 
MF equations and AF energy can be expanded as 



The concave dependence on 8 signals phase separation into a mixed phase of AF and 
PM by the Maxwell construction. The PS extends from the AF phase with density 
n = 1 to a PM phase. Dynamical mean field theory |86] find that for particle or 
hole doping 5 = 1 — n, phase separation occurs between an AF at half filling and 
a paramagnet at a critical filling 8 S depending on the coupling. It has a maximum 
8 S ~ 0.14 around the coupling U ~ 9f where also Tn is maximal. PS is not observed 
in cuprates. It may be inhibited by long range Coulomb repulsion between electrons 
in cuprates as will be discussed later. 

5.2.3 d-wave superfluidity in 2D 

d-wave or d x 2_ Y 2 pairing is particular favoured on a square lattice with on-site re- 
pulsion (which disfavours s-wave pairing) and nearest neighbour attraction as e.g. 
provided by the super-exchange coupling in the Hubbard and Heisenberg model 
with two spin components. The discussion below emphasizes how 2D and also 3D 
optical lattices can resolve important issues as the existence of phase separation, 
stripe phases and dSF 

Repulsive onsite interactions U > inhibit s-wave pairing unless a longer range 
attraction is added such as a nearest neighbour interaction VY,(ij) n i n j- Super- 
exchange in the Hubbard model generates a similar nearest neighbour (spin-spin) 
interaction with coupling J = At 2 /U ', which is believed to be responsible for high 
temperature superconductivity 11871 . In the limit U — > °° extended s-wave pairing 
can occur but requires a strong nearest neighbour attraction V < — n 2 t/2 l74l[88ll . 
whereas d-wave pairing occurs naturally in 2D for even weak nearest neighbour 
attraction. The d-wave mean field gap equation is (see, e.g., IT741I881 ) 




+ 



(83) 



4 V 2E q ' 



(84) 



where now E q = y £q + 4 J rj 2 and r] q = 2 [cos q x — cos q y ] . At low filling the d-wave 
gap can be calculated within mean field to leading orders in the density l90l 
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7tn z 



t 

— +C() + Cl« + C 2 «' 



n< 1, 



(85) 



The higher order corrections in density are: cq = A/n — 1 ~ 0.27, c\ = n/2 — 1 ~ 
0.57 and c 2 ~ 0.09. 

At half filling we can calculate the d-wave gap within mean field as above if 
correlations can be ignored, which requires that the on-site interaction is small. To 
leading logarithmic orders the d-wave pairing gap is 



A d = -yexp -7Zy>7/\V 



n = l. 



(86) 



where again the singular level density leads to a gap similar to Eq. d75l l. It is a 
coincidence that the prefactor 8/e 2 is the same as in Eq. d73l l. Correlations are, 
however, expected to suppress the d-wave mean field gap of Eq. (l86l l near half filling 
when the on-site repulsion is strong. 

T c can in 3D be calculated from the gap equation (|7TT i whereas in 2D a Berezinskii- 
Kosterlitz-Thouless transition occurs ||92l at a lower temperature Tbkt- At low 
density the transition temperature is in 3D generally proportional to the gap with 
the same prefactor T c /A = e Ce j% ~ 0.567 for both s- and d-wave superfluidity 
with and without induced interactions ll28ll . In 2D a mean field critical temperature 
Tmf > Tbkt 11921 . can be calculated from the gap equation as above and we find that 
the same ratio applies even near half filling where logarithmic singularities appear. 
This implies that induced interactions change the mean field critical temperatures 
by the same factor as the gap of Eq. (?). 



5.2.4 t-J-U model predictions 

The various competing phases can be studied in the 2D t-J-U model within the 
Gutz wilier projection method and renormalized mean field theory (RMFT) l90l . 
This method approximates the strong correlations and generally agrees well with 
full variational Monte Carlo calculations [93, 94 1. RMFT predicts phase separation 
near half filling between an antiferromagnetic (AF) Neel order and a d-wave super- 
fluid (dSF) phase for sufficiently strong onsite repulsion. 

The t-J-U model was employed by Laughlin, Zhang and coworkers ||94l to 
study AF, HTc and "gossamer superconductivity" in cuprates and organic super- 
conductors. Both the Hubbard and t-J models are included in the t-J-U Hamiltonian 
H=H v +H t +H s ox 

H = U'£n tt n ii -t £ a] a a ]a + J £ S,Sy , (87) 

where a, a is the usual Fermi creation operator, a = (t,l) is the two hyperfine states 
(e.g. (-§,-!) for 40 K), n ia = a] a a ia the density, S; = "L a & a] a O aa ia i& and {if) 
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Fig. 11 Order parameters for dSF (4sc)> AF (m) and double occupancy (d) vs. doping 5 = 1 — n 
for J ft = At/U = 1/3 in RMFT for the 2D t-J-U model. From (90). 

denotes nearest neighbours. U is the on-site repulsive interaction, t the nearest- 
neighbour hopping parameter and J the spin-spin or super-exchange coupling. 

The t-J-U model allows for doubly occupied sites and thereby also MI transitions. 
As both are observed in optical lattices the t-J-U model is more useful as opposed 
to the t-J model which allows neither. For large U ft the t-J-U and Hubbard models 
reduce to the t-J model with spin-spin coupling J = At/U due to virtual hopping. 
At finite J and U the t-J-U model is to some extent double counting with respect to 
the Hubbard model with / = 0. However, when RMFT is applied the virtual hopping 
and resulting spin-spin interaction vanishes in the Hubbard (t-U) model, i.e. J = 0, 
which justifies the explicit inclusion of the spin Hamiltonian as done in the t-J-U 
model. 

The AF phase dominates near half filling suppressing the superfiuid order param- 
eter. The AF order parameter is m = (vj/2) yjl — T/Tn at temperatures just below 
the critical Nel temperature Tn = 2J. The d-wave superfiuid gap depends sensitively 
on coupling, magnetization and double occupancy as seen in Fig. (JTTJ. It competes 
with the AF order parameter and dominates slightly away from half filling as seen 
in Fig. dm 

However, near half filling phase separation is found between an AF phase at 
half filling coexisting with a dSC phase at a density |x|^S0.12 as determined by the 
Maxwell construction. The PS is found to terminate at coupling J ~ 0.55f (U ~ 7.3f ) 
where the double occupancy undergoes a first order transion from zero to a finite 
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Fig. 12 Density distributions of Fermi atoms in a zero temperature 2D optical lattice confined by 
a harmonic trap for U = 12/ (7 = r/3) filled with the number of particles N/N c = 0.5, 1, 5, and 10. 
The density discontinuities from n ~ 1.14 and n ~ 0.86 to the antiferromagnetic insulator at n = 1 
due to PS between dSF and AF are absent for stripe phases (see text). Thus the density distribution 
changes to a Capitol form rather than the usual wedding cake form for Mott insulators. 

value. The PS has a curious effect on the density distribution in a trap as shown in 
Fig. ([T2l i. The trap size is of order R c = -J 8f JhG)a 0SC , and the number of trapped 
particles N = 2k J n(r)rdr is of order N c = 7tR^ in 2D. 

The density distributions and MI plateaus have been measured experimentally 
for Bose atoms in optical lattices by, e.g., differentiating between singly and doubly 
occupied sites [71 ). Recently Schneider et al. ||89l have measured column densities 
for Fermi atoms in optical lattices and find evidence for incompressible Mott and 
band insulator phases. Even lower temperatures are required for observing the AF 
and dSF phases and the density discontinuities due to PS. 

Spin and charge density waves in form of stripes are not included in the above 
RMFT calculations. Stripes are observed in several cuprates II 1001 whereas numeri- 
cal calculations are model dependent. Long range Coulomb frustration can explain 
why PS is replaced by stripes and an AF phase at very low doping l90l as observed 
in cuprates. 



5.3 High temperature superconductivity vs. optical lattices 

After Kamerlingh Onnes discovered superconductivity in 1911 it took nearly half a 
century before Bardeen, Cooper and Schriefer (BCS) in 1957 understood the physics 
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behind the manybody pairing mechanism. In comparison half the time has passed 
since Bednorz & Miiller in 1986 discovered high temperature superconductivity 
(HTc). Although important progress has been made on the mechanisms behind this 
important phenomenon, a full understanding is still lacking. Most agree that HTc is 
d-wave and can probably be described by the one band Hubbard model on a two- 
dimensional (2D) lattice with strong repulsive on-site interaction U ~ 3f, but no 
consensus has been reached on the origins and influence of the pseudogap, charge 
and spin gaps, stripes, etc. 

A sketch of a typical phase diagrams of HTc cuprates is shown in Fig. (fT3l l. 
The antiferromagnetic phase extends around half filling and is asymmetric due to 
next neighbour hopping terms t '. The Neel temperature drops rapidly from its max- 
imal value 7)v(5 = 0) - somewhat lower than the Fermi temperature. The maximal 
J C ^,100K are collected in the "Uemura" plot l99l . It is typically between one and 
two orders of magnitude lower than the Fermi energy which again is several orders 
of magnitude larger than conventional superconductors as Nb, Sn, Al, Zn. 

Optical lattices hold great promises for solving the HTc problem [95, .96] since 
densities, temperatures, interactions, spins, etc. can be varied in a controlled man- 
ner and studied in experiments, where the sub nanometer lattice constants in solids 
are upscaled to micrometers in optical lattices. One serious experimental hindrance 
remains, namely the extremely low temperatures required, but we have hopes that 
the experimentalists can overcome this problem rapidly as in so many other cases 
with cold atoms - and it certainly does not hinder theoretical studies. 

5.3.1 Stripes 

Stripes have been discovered in low-energy magnetic neutron scattering in doped 
cuprates at incommensurate longitudinal and horizontal charge and spin wave num- 
bers COQl, e.g. Q £ . = (27r/A)(0,±<5/<S.v) and Q s = (n/k)(l, 1 ± S/S s ) for the lon- 
gitudinal charge and spin wave numbers respectively. Here, the average doping or 
holes are all in the stripes with hole filling 8 S embedded in an AF with filling no holes 
(half filled) such that the average hole density is 8 = 1 — n. The stripes are half filled 
8 S = 1/2 for low average dopings 8 < 1/8 but at larger doping it increases linearly 
to filled stripes 8 S = 1 such that 8/8 s = 1/4 remains constant. Therefore the charge 
(spin) density stripes appear with periodic distance d {2d) depending on doping as 
d = X8 s /8. The stripe distance decreases with increasing doping until 8 > 1/8, 
where after the stripes remain at a distance d = 4 A . The stripes act as anti-phase 
domain walls and the spin density wave therefore has periodicity twice the length 
of the charge-density wave. When 1/8 < 8 < 1/4 the 2D spin orientation therefore 
looks like 

...tit°m°tit°m°- 
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Fig. 13 Typical phase diagram of high temperature superconductors. Neel temperatures (7]v) and 
critical temperatures for d-wave superconductivity (T c ) are plotted vs. particle (5 < 0) and hole 
(5 > 0) doping. T* is the pseudogap (see text). 



etc., in case of vertical stripes o. They are spin-balanced at filling 1 — S s and occur 
with periodicity d = AX . 

Diagonal, horizontal/vertical, checkerboard stripe solutions have been found in a 
number of models. In MF models Zaanen and Gunnarson l97l found stripes of hole 
density S s = 1 that are vertical or horizontal for U /f^S3 —4 and diagonal otherwise. 
In DMRG calculations l98l stripes with hole density S s = 1/2 are found in agree- 
ment with experiments 11001 . Stripes appear at temperatures above HTc but below 
the pseudogap T* 11011 and there are some evidence that charge density waves ap- 
pear at slightly higher temperatures than spin density waves. 

The stripe phase is a specific ordered mixed AF and PM phase and is a continuous 
transition between the two pure phases as function of density. Similar mixed phase 
solutions are believed to occur in neutron star crusts between nuclear matter and 
a neutron gas and possibly also between quark and nuclear matter 11021 . In both 
cases Coulomb energies add complexity to the mixed phases by ordering them into 
structured crystallic phases. 

The various MI, AF, stripe, d- and s-wave superfluid phases have distinct density 
and momentum correlation functions |fT^ rT05l [7T1 [72l [731 rT06l [90ll . 
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5.3.2 Coulomb frustration 



Long range Coulomb interactions prevent phase separation into two bulk phases 
of different charge density but may not inhibit the formation of localized holes, 
pairs and stripes. The phase diagrams of cuprates with such Coulomb frustration 
and optical lattices without may therefore be very different. Coulomb frustration in 
cuprates has been discussed in connection with stripes (see e.g. 11103111041 ). 

In the following we shall consider the stripes as rods with charge less than the 
surrounding phase and calculate the additional Coulomb energy of such structures. 
Coulomb energy densities have been calculated for structures of various dimension- 
ality D and volume filling fraction / 110211 



E C = —Ap 2 R 2 f 
D + 2 F ' 



(i-^-^+z 



D-2 X 2 



(88) 



Here, the charge density difference between the two phases is Ap = e5 s /A 3 for the 
stripes and the volume filling fraction is / = 8/8 s = X/d. The dimensionality is 
D = 3 for spherical droplets or bubbles, D = 2 for rods and tubes and D = 1 for 
plate-like structures. The diameter of the spheres, rods or the thickness of the plates 
is of order the distance between layers 2R ~ A ~ 4 A. The stripes are rods in a 2D 
plane but are embedded in a 3D layered structure. For D = 2 the expression in the 
square bracket of Eq. ( f88l > reduces to [ln(l//) — 1 +/]. The logarithm originates 
from the Coulomb integral / dz/z along the rod length z, which is cutoff by other 
rods at a length scale I ~ A s/J. For the cuprates we furthermore reduce the Coulomb 
field by a dielectric constant of order Ed ~ 5. The resulting Coulomb energy of stripe 
or rod-like structures D = 2 is 

ec^f^/[ln(l//)-l+/]. (89) 

Energy costs associated with the interface structures are usually added. Such 
surface energies are difficult to calculate for the stripes because their extent is only a 
single lattice constant. In principle they are already included in the stripe models. We 
will therefore just add the Coulomb energies given above. However, the Coulomb 
energies and the energy of the systems as a whole, may be reduced by screening and 
hole hopping into the AF whereby R increases but Ap and / are reduced. 

Inserting numbers e 2 /he = 1 /137, Ed = 5, 8 S = 1/2, A = 4A, we find that E c ~ 
150 me V ( 5/A 3 ) [In ( 1 //) — 1 + /] . In comparison the energy gain by changing phase 
from an AF to a stripe or PM phase increases with doping as ~ /5/A 3 , where 
J ~ (0. 1 - 0.2)? ~ 50 - 100 meV in the cuprates. The Coulomb energy of Eq. (89} 
thus dominates at small doping due to the logarithmic singularity and therefore the 
AF phase of density n^l is preferred. 

The AF phase is the ground state as long as the Coulomb energy of Eq. 
exceeds /a/A 3 corresponding to doping less than 
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Inserting the above numbers we find 8af — 0.1 which is within range of the observed 
\8af \ — 0.03 for hole doped and 8af — 0.15 for particle doped cuprates. In MF the 
particle-hole asymmetry arises from the next-nearest neighbour hopping t' ~ — 0.3? 
and leads to an AF phase extending from half filling up to a particle doped density 

n> 1. m 

In the above picture the incommensurate stripe phases at small doping arise due 
to Coulomb frustration when t' = 0. At larger doping the stripes approach each other 
and will eventually affect each other. Experimentally the stripes undergo a transition 
from an incommensurate to a commensurate phase at 8 ~ 1/8 corresponding to a 
stripe periodicity of four lattice spacings. 

The Mermin- Wagner theorem states that a continuous symmetry cannot be spon- 
taneously broken at finite temperature in two and lower dimensional systems for suf- 
ficiently short range interactions |9"T1l , and as a consequence there can be no phase 
transition at finite temperature. In stead a Berezinskii-Kosterlitz-Thouless transition 
can occur as found in the 2D Hubbard model (92). This is incompatible with the 
phase diagram of HTc cuprates where AF, dSC and stripe phase transitions are ob- 
served. These paradoxial transitions and the origin of the phases are still unresolved 
problems which are largely avoided in the literature. Some suggest that inter-planar 
couplings effectively increase the dimension above two. Others that correlations 
lengths are so long that they effectively looks like real transitions. Another sugges- 
tion in line with the Coulomb frustration discussed above is that Coulomb interac- 
tions are responsible for the phase transitions to stripes and d-wave superconductiv- 
ity. The long range Coulomb interactions are not comprised in the Mermin- Wagner 
theorem, which only applies to sufficiently short range interactions U(r) °< r~ a with 
a >D + 2 (9D. 



6 Summary and Outlook 

In the last decade the physics of cold atoms has brought important understanding 
of the unitarity limit, universality, the BCS-BEC crossover, and many other prop- 
erties of cold Fermi atomic systems as describe above and elsewhere in this vol- 
ume. Presently a number of other interesting phenomena are investigated such as 
the three component system of 6 Li and other multicomponent systems, the ferro- 
magnetic transition at large positive (repulsive) scattering lengths, spin polarized 
systems, phases in optical lattices, etc. 

In these lecture notes the main topic has been universality in crossover. Not only 
the BCS-BEC crossover in uniform system but also in the repulsive "ferromagnetic" 
crossover, in multicomponent systems, traps and lattices. The Fermi particles have 
mostly been atoms but applications to neutron, nuclear and quark matter, nuclei and 
electrons in solids have been made wherever possible. 
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In view on the important impact that the question Bertsch posed a decade ago 
made, we may attempt to ask new relevant questions: 

What are the universal thermodynamic functions and parameters for positive 
scattering lengths and is there a ferromagnetic phase transition as a — > +°°? 

- in both uniform systems and traps. One opinion is given above but a dispute is 
ongoing [12] [H. 

A most promising direction is optical lattices where we may ask: 

What are the phase diagrams of the two and three dimensional Hubbard mod- 
els as realized in optical lattices? 

Attractive interactions will lead to a superfluid state with critical temperatures that 
have been calculated in detail by Monte Carlo. For repulsive interactions the phase 
diagram is not very well known in 2D or 3D even for the single band t-U or t-J 
model. The antiferromagnetic phase may extend to densities near half filling before 
undergoing a transition to a paramagnet but it is competing with phase separation, 
stripe phases and a ferromagnetic phase for very strong repulsion. In 2D it also 
competes with d-wave superconductivity. 

A related question, which is particular relevant for HTc, is: 

Does the 2D Hubbard model (as realized in optical lattices) exhibit high tem- 
perature superconductivity (as seen in cuprates)? 

Because the electrons in HTc have charge they cannot undergo phase separation in 
bulk, which in turn may be responsible for stripe formation in cuprates in the region 
of densities and temperatures where also d-wave superconductivity appears. Ultra- 
cold Fermi atoms are neutral and can therefore undergo phase separation near half 
filling as is predicted in some models, e.g. the t-J-U model with Gutzwiller projec- 
tions discussed above. If so, the region of densities near half filling with antifer- 
romagnetic and d-wave superconductivity do not exist in optical lattices. Cuprates 
may also have multiband coupling and certainly has next nearest neighbour hop- 
pings, which cause asymmetries in the phase diagrams around half filling for both 
the AF and dSC phases. In any case, the differences between the phase diagram 
of the cuprates and the Hubbard model as realized in optical lattices will reveal 
the important driving mechanisms and possible additional ingredients present in the 
phenomenon of HTc. 2D optical lattices may also resolve the Mermin- Wagner para- 
dox. 

In summary, we have good reasons to believe that ultracold Fermi atoms in traps 
and optical lattices will in few years bring deeper understanding of strongly corre- 
lated systems, its para-, ferro- and antiferromagnetic phases, Mott insulators, stripes, 
s- and d-wave superfluidity, and high temperature superconductivity. 
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